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Understanding heat generation and transport processes in a magnetic tunnel junction (MTJ) is 
a significant step towards improving its application in current memory devices. Recent work has 
experimentally demonstrated the magneto-Seebeck effect in MTJs, where the Seebeck coefficient of 
the junction varies as the magnetic configuration changes from a parallel (P) to an anti-parallel (AP) 
configuration. Here we report the study on its as-yet-unexplored reciprocal effect, the magneto- 
Peltier effect, where the heat flow carried by the tunneling electrons is altered by changing the 
magnetic configuration of the MTJ. The magneto-Peltier signal that reflects the change in the 
temperature difference across the junction between the P and AP configurations scales linearly 
with the applied current in the small bias but is greatly enhanced in the large bias regime, due to 
higher-order Joule heating mechanisms. By carefully extracting the linear response which reflects 
the magneto-Peltier effect, and comparing it with the magneto-Seebeck measurements performed 
on the same device, we observe results consistent with Onsager reciprocity. We estimate a magneto- 
Peltier coefficient of 13.4 mV in the linear regime using a three-dimensional thermoelectric model. 

Our result opens up the possibility of programmable thermoelectric devices based on the Peltier 
effect in MTJs. 


The electrical resistance of a magnetic tunnel junction 
(MTJ), a stack of two ferromagnetic layers separated by 
an insulating tunnel barrier, depends on the relative mag¬ 
netic orientation of the two magnetic layers [1-3]. This 
tunnel magnetoresistance (TMR) effect puts MTJs at the 
forefront of the applications in the field of spintronics [4] . 
Spin caloritronics [5-7] is an emerging field that couples 
thermoelectric effects with spintronics. Many interest¬ 
ing physical phenomena were discovered such as the spin 
(-dependent) Seebeck effect in magnetic metal [8], mag¬ 
netic semiconductor [9] and magnetic insulator [10]. Par¬ 
ticularly, in spin tunneling devices, the magneto-Seebeck 
effect was theoretically studied [11-14] and experimen¬ 
tally observed [15-19] in MTJs, where the Seebeck coef¬ 
ficient of the junction can be varied by changing the mag¬ 
netic configuration. More recently, the spin (-dependent) 
Peltier effect that is driven by spin (polarized) currents 
has been experimentally observed in metallic [20, 21] and 
insulating ferromagnets [22] , which are shown to obey the 
Thomson-Onsager reciprocity relation [23-25] to the spin 
(-dependent) Seebeck effect. From this relation, the re¬ 
ciprocal effect of the magneto-Seebeck effect, which can 
be named as magneto-Peltier effect, is also expected in 
MTJs (see Fig. l(a)(b)). 

However, experimental studies of the magneto-Peltier 
effect have not been reported so far. Its small effect com¬ 
pared to the often-dominant Joule heating effects has left 
the experimental observation elusive. In this work, we re¬ 
port the first experimental study of the magneto-Peltier 
effect as well as higher order heating effects, and compare 
the Peltier measurements to the Seebeck measurements 
on the same junction. Via sensitive thermometry archi¬ 
tecture and measurement techniques, we are able to mea¬ 


sure small temperature changes as well as distinguish the 
linear (due to Peltier effect) and nonlinear effects (due to 
Joule heating). 

Although both the electric conductance and Seebeck 
coefficient depend on the relative magnetic configura¬ 
tion in a MTJ, the mechanisms behind them are not the 
same. While the electric conductance is determined by 
the transmission probability Tp^ap{E) of electrons across 
the insulating barrier around the Fermi energy Ep, the 
Seebeck coefficient Sp^ap solely depends on the electron- 
hole asymmetry of Tp^ap{E) around Ep. By Onsager 
reciprocity, the Peltier coefficient H is closely related to 
S' by n = STq, where Tg denotes certain temperature. 
Using the expression for S [11, 15] we can express H as: 

n jTp,Ap{E)iE-Ep){-dEfo)dE 

ejTp^Ap{E){-dEfo)dE ’ ^ ^ 

where e is the elementary charge and /g is the Fermi- 
Dirac distribution function at temperature Tg. For a con¬ 
stant tunneling current / through the MTJ, the Peltier 
heat Qu carried by this current is different for the P 
and AP configurations, leading to different temperature 
biases across the junction between the two configura¬ 
tions. Disregarding Joule heating and assuming that 
the thermal conductance k of the MTJ is independent 
of the magnetic conhguration, the temperature differ¬ 
ence between the parallel and antiparallel conhguration 
AT = ATap — ATp can be estimated by balancing the 
Peltier heat Qu with the backhow of the heat current 
through the junction as 

AT= ^(H^p-Hp), (2) 
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FIG. 1. Concept and device geometry, (a), Concept of the 
magneto-Seebeck effect. A temperature gradient is applied 
across an MTJ, resulting in a Seebeck voltage that is de¬ 
pendent on the magnetic configuration, (b), Concept of the 
magneto-Peltier effect. A charge current is sent through an 
MTJ, resulting in a Peltier heating/cooling at the interfaces 
that also depends on the magnetic configuration. Joule heat¬ 
ing is not shown for simplicity, (c), Schematic representation 
of the measured device. A Pt-NiCu (constantan) thermo¬ 
couple that is electrically isolated from the top contact by 
an AI 2 O 3 layer (green color) is used to detect temperature 
changes. In the Peltier measurement, charge current is sent 
through the pillar (from contact 1 to 2 ), while recording the 
voltage over the thermocouple (contact 3 and 4), as plotted 
here. Contacts 5 and 6 are used for 4-probe TMR measure¬ 
ments. For the reciprocal Seebeck measurement, current is 
sent through the thermocouple while recording the voltage 
over the pillar using contacts 1 and 2. (d). Optical microscope 
image of the measured device. The dotted circle indicates the 
location of the MTJ. The size of the junction measured in the 
main text is 2.7x5 fim^ by size. 

where t is the thickness of the tunnel barrier and A is the 
area of the junction. Using the parameters given in ref. 
[15], it can be estimated that for an electric current den¬ 
sity of 5x10^ A/cm^ through the junction, the change in 
temperature due to the magneto-Peltier effect can reach 
^100 fjK at room temperature, which requires sensitive 
thermometry technique, as we use here. 

Fig. 1(c) (d) show the device geometry and measure¬ 
ment conhguration as employed in our experiment. We 
study the CoFeB/MgO/CoFeB MTJs, same as in ref. [15] 
that are reported to have a large magneto-Seebeck effect 
[11, 15, 16] due to their half-metallic transmission prop¬ 
erty [3]. The layer structure of the patterned MTJ sample 
stack (from bottom to top) is: Si ( 100 )/SiO 2 500 nm/Ta 
15 nm/CoFeB 2.5 nm/MgO 2.1 nm/Co-Fe-B 5.4 nm/Ta 
5 nm/Ru 3 nm/Cr 5 nm/Au 25 nm. Detailed fabrica¬ 
tion processes can be found in ref. [15]. To sense the 
temperature change locally at the top of the junction, a 
thermocouple consisting of constantan (Ni 45 Cu 55 ) and Pt 
is fabricated over the top contact of the MTJ [21, 24, 26]. 
Constantan is dc sputtered and Pt is e-beam evaporated. 


both 90 nm in thickness. A 50-nm-thick AI 2 O 3 layer is 
e-beam evaporated over the top contact to electrically 
isolate MTJ and thermocouple from the top contact, so 
that no charge-related effects are picked up by the ther¬ 
mocouple. Finally, a 130 nm layer of Au is deposited to 
connect the two arms of thermocouple, creating a uni¬ 
form temperature distribution over the junction. In the 
Peltier measurement configuration, we send a charge cur¬ 
rent through the pillar (contact 1 to 2 ) while record¬ 
ing the thermovoltage using the thermocouple (contact 
3 to 4). Note that in this configuration the temperature 
changes resulting from both the Peltier and Joule heat¬ 
ing effects are measured. Meanwhile we also record the 
tunnel magnetoresistance (TMR) by four-probe method 
using contacts 5 and 6 . All measurements shown in the 
main text are performed on a single device at room tem¬ 
perature. Measurements on two other samples can be 
found in the supplementary. 

A standard lock-in technique is used for our measure¬ 
ments, where an ac current is sent through the system at 
a low excitation frequency (~17 Hz), so that a steady- 
state temperature condition is reached and at the same 
time capacitive coupling is suppressed. The voltage out¬ 
put from the sample is separated into different harmonic 
signals in terms of the input frequency. 

With this technique it is possible to isolate the Peltier 
effect that is linear with current from the Joule heating 
effect that is of second or even higher responses. In a 
simple case where only and are present, Peltier 
effect is linked to while the secondary Joule heating 
effect is linked to . But in a more complex case where 
higher order effects are also present, the n-th order signal 
is not directly equal to the n-th harmonic signal; instead, 
all higher harmonic signals (with the same parity as n) 
that are nonzero need to be included by a straightforward 
algebraic operation (see supplementary for more details). 

Fig. 2 shows the experimental results of the magneto- 
Peltier measurements. We apply an ac current of 150 
/iA (r.m.s) through the MTJ while sweeping the mag¬ 
netic field. Both the first and second-harmonic voltages 
recorded at the thermocouple (U^^ and shown in 

Fig. 2(a)(b), exhibit four abrupt changes corresponding 
to the switching from P to AP configuration and back, 
implying a change in the temperature at the top con¬ 
tact. of TMR is shown in Fig. 2(c). We measure the 
current dependence of these signals correspondingly, for 
both the P and AP configurations, as shown in Fig. 2(d)- 
(f). We did not apply ac currents higher than 150 fiA 
to avoid dielectric breakdown of the MTJ (~2 V across 
the junction for 2.1 nm MgO) [27]. In the P configu¬ 
ration we have a simple case where no higher harmonic 
signals can be detected, and the signal de¬ 

tected at the thermocouple is linear (quadratic) with the 
current, which can be considered as the Peltier signal 
(Joule heating signal). However, for AP configuration. 
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FIG. 2. Lock-in signals from magneto-Peltier and TMR 
measurement configurations. (a)(b), First {V^^) and sec¬ 
ond (V^^) harmonic signal measured at the thermocouple 
for an r.m.s current of 150 fiA. c, First harmonic sig¬ 

nal of TMR measured at the same current. On the right 
axes the temperature differences detected by the thermo¬ 
couple relative to the room temperature To are given as 
Trc -To = - Smcu), where To is 290 K. (d)- 

(f), Current dependence of the corresponding measurements, 
for P (open triangles) and AP (open circles) configurations, 
where P (AP) is obtained when setting the magnetic field at 
10 (2) mT. Solid lines are linear (in d, f) or quadratic fits (in 
e), as references. 


we have a more complex case where higher order effects 
are present. The I — V characteristic of the MTJ (see 
supplementary) is nonlinear with the current, or in other 
words, the resistance of the junction R is bias-dependent 
[28-30]. Therefore, the Joule heating effect (/ • V) is not 
only present in the second order, but also brings on higher 
order responses. The consequence of this nonlinearity is 
twofold: first, R decreases with both larger positive and 
negative biases. This leads to an even higher order re¬ 
sponses at the thermocouple which deviate the from 
a quadratic behavior, as shown in Fig. 2(e). In addition, 
R also shows an asymmetric dependence for +I and — 
i.e., V{+I) ^ V{—I), indicating that the dissipation at 
the junction is different when the bias is reversed. The 
reason for this asymmetry can be attributed to the in¬ 
evitable difference between the two interfaces across the 
MgO [30]. Although this effect is only present in higher 
odd order heating signals on thermocouple, it mimics a 
Peltier-like effect and strongly deviates from a linear 


behavior (see Fig. 2(d)). A more extensive quantitative 
analysis can be found in the supplementary. 

It is therefore important to determine the pure lin¬ 
ear signal (i?i/o) of the AP case in order to discuss the 
magneto-Peltier effect. Taking advantage of the lock- 
in detection technique, we can measure the higher har¬ 
monic signals ...) by tuning the lock-in detec¬ 

tion frequency to the corresponding harmonics frequency. 
The first order (linear) response can be expressed as 
a linear combination of the odd harmonics voltages as 
RiIq = + 3V^^ + 5V^-I'... (see supplementary). Here 

we only include higher harmonic signals up to as 

cannot be determined accurately within our noise 
level (< |5 nVj). The results are shown in Fig. 3(a). 
However, the difference between P and AP still shows 
a seventh order behavior as a function of the current. 
This means that, especially at larger currents (above 80 
fiA), the seventh order response is still present; however, 
the P^^ signal that can be measured is only 1/8 of the 
seventh order response, which made it difficult to be in¬ 
cluded to extract RiIq- Nevertheless we can still rely on 
the lower-current regime before the onset of the seventh 
order response (below 80 fj,A, circled part in Fig. 3(a)), 
which can be regarded as purely linear regime. We fit the 
curves for P and AP individually for this regime and es¬ 
pecially focus on their difference, which can be considered 
as the magneto-Peltier effect. Although the difference is 
small compared to the noise, we fit it linearly and esti¬ 
mate a slope range bounded by one standard deviation, 
58±35 as shown in the inset of Fig. 3(b). 

To support our estimation of the magneto-Peltier sig¬ 
nal, we perform magneto-Seebeck measurements on the 
same device by effectively reversing the role of the current 
and voltage contacts as used in the Peltier measurement. 
Here we send an ac current through the thermocouple 
(contact 3 to 4) thereby creating a vertical temperature 
gradient over the MTJ via the Peltier heating/cooling 
(oc I) at the NiCu-Au and Au-Pt interfaces. The Seebeck 
voltage (open-circuit thermovoltage) due to this vertical 
temperature gradient is measured using contacts 1 to 2. 
In Seebeck measurements, it is possible to send larger 
currents up to 2 mA through the thermocouple with a re¬ 
sistance of 190 H. Unlike the Peltier measurement, in the 
Seebeck measurement no higher odd harmonic features 
are observed for either P or AP configurations, implying 
a linear behavior for the Seebeck signal in the measured 
current range. This is because the thermocouple is purely 
ohmic, in contrast to the nonlinear MTJ. The current 
dependence of the magneto-Seebeck measurements are 
shown in Fig. 3(c). According to the Thomson-Onsager 
reciprocity relation, the linear response signals for the 
Peltier and Seebeck effect should be the same, as well 
as the difference between P and AP configurations [24]. 
From Fig. 3(c), the difference between the two config¬ 
urations is 12.5±0.4 /rfl, which falls into the estimated 
range of the magneto-Peltier effect within 2 standard de- 
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FIG. 3. Comparison of the magneto-Peltier and magneto-Seebeck linear responses, (a), The extracted linear response calculated 
up to including the , as a function of applied current. The circled part is where the is close to 0 and is therefore represents 
pure linear signals. Inset shows the difference between the P and AP configurations which is htted to a seventh order behavior, 
(b), Linear fitting of the circled part in (a). The shaded zones indicate the standard deviations of the fitted slopes. Inset 
shows the difference between the P and AP configurations which is linearly fitted, with the shaded area indicating one standard 
deviation, (c), Current dependence of the Seebeck measurements, without any compensations from higher harmonic signals. 
Inset shows the difference between the AP and P configurations which is htted to a linear behavior. 


viations (corresponds to a confidence level of 95%), there¬ 
fore showing no statistically significant difference. This 
is consistent with the reciprocity between the magneto- 
Seebeck and magneto-Peltier measurements. In our opin¬ 
ion, there is no fundamental reason for the rather large 
difference in the average values for magneto-Peltier and 
Seebeck coefficients, except for the experimental difficul¬ 
ties in obtaining the magneto-Peltier coefficient. Note 
that the background signals for Seebeck and Peltier con¬ 
figurations correspond closely, indicating the validity of 
our approach. However, the backgrounds contain See- 
beck/Peltier effects from all metal interfaces, and there¬ 
fore are not directly linked to the Seebeck/Peltier coeffi¬ 
cients of the MTJ. 

By using a three-dimensional finite element model (3D- 
FEM) [31] we can quantify our results. We focus on the 
estimation of the relative change of the Peltier coeffi¬ 
cient from P to AP configuration of the MTJ. We do not 
model the electron tunneling process, but regard MgO 
as a conductor whose electrical conductivity and Peltier 
coefficient vary between P and AP states, while keeping 
other properties of the MTJ constant. The details can 
be found in supplementary. We find that the modeled 
magneto-Peltier signal is very sensitive to the choice of 
the thermal conductivity of the MgO layer (KMgo) and 
the difficulty of measuring this quantity directly can cre¬ 
ate a big uncertainty in our estimation. Here we adopt 
the same value from ref. [15], where KMgO = 4 W/(m-K) 
was used for 2.1 nm MgO layer, taking into account 
both the crystalline quality of MgO and its thermal in¬ 
terfaces with CoFeB. By fitting to our experimental re¬ 
sult 12.5±0.4 from the Seebeck measurements (which 
has less statistical uncertainty), we obtain the change of 
Peltier coefficient of MgO to be AH = H^p—Hp = AS-T 


=13.4 mV from P to AP, which by Onsager relation cor¬ 
responds to AS = Sap — Sp =46.2 /rV/K. This is close 
to the Seebeck coefficient change of MgO reported in ref. 
[15]. 

In conclusion, we have observed the magneto-Peltier 
effect in magnetic tunnel junctions and confirmed its reci¬ 
procity to the magneto-Seebeck effect by measuring both 
effects in a single device. We also observed higher or¬ 
der heating effects which greatly enhance the magneto- 
Peltier signal in the large-bias regime. We attribute this 
effect to the asymmetric resistance of MTJ for opposite 
bias. In addition to providing additional insight in the 
nature of heat dissipation in MTJs, our results open up 
the possibility of magnetically controllable cooling mech¬ 
anism in MTJs, which can be potentially applied in novel 
magnetic logic devices. 
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SUPPLEMENTAL MATERIAL 

I. Results for two other samples 

In the main text, all the measurements shown were performed on a single device, which we here denote as sample 
A. Here we show additional measurements performed on two other samples, sample B (Fig. 4(a)-(c)) and sample C 
(Fig. 4(d)-(f)). From left to right the first harmonic , second harmonic signal of magneto-Peltier measurements 
and of TMR measurement are plotted. Similar features are observed compared to Fig. 2 in the main text. 


_ sample.B_ 



magnetic field (ml) 



(c) 



magnetic field (ml) 



FIG. 4. Magneto-Peltier and TMR measurements for two other samples. Figure (a)(b)(c) are measured on a sample where the 
junction size is 2.1x3.2 pm^. Figure (d)(e)(f) are measured on a sample where the junction size is 1.8x3.1/rm^. Here the r.m.s 
current for both measurements is 50 pA. 


II. Extraction of the linear response from the lock-in measurements 

The measurements shown in this paper are all performed using a lock-in detection technique. In this technique, an 
ac current I{t) = V^Io sin(a;t) with an angular frequency w, generated by one of the lock-in amplifiers, is sent to the 
studied system as an input current. In general, the system output in terms of voltage can be written as: 

V(t) = i?i/(t) -f i?24^(t) -|- i?3/^(t) -|- -|- i?5/^(t) -|- ... (3) 

where is the n-th order response and is the n-th order coefficient, which we want to determine. The 

lock-in technique can separate different harmonic contributions of V (t) according to the orthogonality of sinusoidal 
functions: it multiplies V(t) by a reference sine-wave signal, with an angular frequency that is an integral (denoted 
by n) multiple of w, and then averages this product over time as: 

^ f sin{nu}t + (j)))V{t)dt, (4) 

where we can obtain a signal as the n-th harmonic signal directly from the output of the n-th lock-in amplifiers. 
In this way, different harmonic components of V(t) are separated. However, the n-th harmonic signal is often 
not equal to the n-th order response i?„/"(t), especially at large current biases where higher order responses are 
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FIG. 5. Higher harmonic signals for magneto-Peltier measurement. Magnetic field sweep with Peltier measnrement configuration 
for an r.m.s current of 150 /rA, where (a)I/^^, (c)H®-^ and {d)V’^^ are shown. 


non-negligible, although it is possible to express one from the other. Assuming we have a voltage response only up to 
the fifth order, different harmonic signals can be expressed as: 


= RJo + (0 = 0“), 

(5a) 

y2/ = ^(i?2/o" + 2i?44) (</> = -90°), 

(5b) 

p3/ = _l(2i?3/3 + 5i?5/o') (<^ = 0“), 

(5c) 

V^f = -^RiI^ (<^ = -90°), 

(5d) 

V^f=^-R,4 (0 = 0°). 

(5e) 


It is then possible to find the n-th order response as a linear combination of , for example, for linear response, 
we have 

RJo = . ( 6 ) 

In the magneto-Peltier measurement, we are interested in RiIq, which reflects the (magneto-)Peltier effect. By 
using three lock-in amplifiers to record and simultaneously, it is possible to determine RiIq using 

Eq. ( 6 ). Such calculation has been performed to obtain Fig. 3(a) of the main text. 

Fig. 5 shows different higher harmonic signals and V'^-^, measured with an r.m.s current of 150 fiA 

on the same sample as main text, where the higher harmonic signals are the most prominent in our measured current 
range. For the odd harmonic signals, switches from P to AP are visible in and , but not in V^f, as the signal 
is already within the noise level. Therefore, it is reasonable to perform Eq. ( 6 ) to obtain linear response in our results. 

For the even harmonic signals, they are only relevant if we want to obtain even order responses. Here we present 
the magnetic field dependence of as shown in Fig. 5(b). Taking the contribution of into consideration by 
using Eq. (5b), we can get the second order response R 2 I 0 whose dependence on the applied current is close to a 
quadratic behavior. 





















III. Finite element thermoelectric model 


To understand the thermoelectric transport properties and extract the Peltier coefficient of the MTJ we make use 
of a finite element model [31] that solves the steady-state three-dimensional current and heat equation. In this model 
the charge current density J, driven by the applied voltage, and the heat current Q, driven by the thermal gradient, 
are related to each other as 



f a aS\ 
l^crn K ) [yTj 


( 7 ) 


where a is the electrical conductivity, S is the Seebeck coefficient, 11 = STq is the Peltier coefficient at any temperature 
Tq. Conservation of charge and Joule heating are included as 


V • 



(jVct) (^a{V^V + S'^V^T + 2SVV-VT)) 


( 8 ) 


The material parameters used in the FEM are shown in Table 1. We do not model the tunneling process directly 
but consider the MgO as a material whose electrical or thermal parameters, such as the conductivity and Peltier 
coefficient, depend on the magnetic configuration of the MTJ, a similar approach as employed in Ref. [15, 18, 32]. 
From the four-probe resistance of the MTJ (TMR measurement), we obtain the conductivity for the parallel up and 
antiparallel configurations (tap and assign these conductivity values to the MgO barrier. The thermal conductivity 
of the insulating MgO barrier is estimated to be 4 W/(m-K), according to Walter et al. [15]. The Peltier coefficient 
of MgO is then varied to fit our experimental results. 

In the following, we describe our approach for extracting the Seebeck (Peltier) coefficient of the MTJ. We first 
start by obtaining the electrical conductivity of the MgO from the tunnel magnetoresistance measurements. In device 
A where we observed a TMR of 87.5% (with parallel resistance i?p=1.6 kO and antiparallel resistance i?AP=3.0 
kU) when I = 150 /iA, we obtain the corresponding conductivities crp=0.093 S/m and aAP=0-05 S/m, respectively. 
Although the thermal conductance of MgO is expected to depend on the magnetic configuration of the MTJ [13], for 
simplicity, we do not take this scenario into account. In the estimation of the magneto-Peltier coefficient, we match 
the measured Peltier signals by varying the magneto-Seebeck coefficient AS = ASap — ASp. Using the results from 
earlier studies of Walter et al. in similar MTJs (Seebeck coefficient values Sp= -108 ^V/K and Sap= -99 /rV/K), we 
obtain a Peltier signal AU = AVap — AVp of 0.37 nV for I = 150 ^A, which is 5 times smaller than our experimental 
result (1.9 nV, corresponding to 12.5 /rU). Increasing AS' by 5 times to 46.2 /rV/K fits our result and gives a AH of 
13.4 mV for the magneto-Peltier coefficient of the MTJ. 


IV. Joule heating in higher order responses 

In the magneto-Peltier measurements as shown above, non-negligible higher odd harmonic signals have been ob¬ 
served, as plotted in Fig. 5. In this section, we show that these signals can be ascribed to Joule heating related 
mechanisms. 


TABLE I. For Au, Pt and NiCu, a was measured in dedicated devices [33] and for Ta, a is obtained from Ref. [34]. The 
thermal conductivities k were obtained using the Wiedemann-Franz relation valid at the temperature of our experiment. The 
rest of the material parameters were taken from Walter et al. [15]. For MgO, cr is determined from TMR measurement and S 
is varied as a fitting parameter. 


Material 

(thickness) 

cr 

( 10 ® s/m) 

S 

(hV/K) 

K 

(W/(m.K)) 

Au (130 nm) 

27 

1.7 

180 

Pt (90 nm) 

4.8 

-5 

37 

NUsCuss (90 nm) 

2 

-30 

20 

AI 2 O 3 (50 nm) 

0 

- 

0.15 

Au top contact (25 nm) 

18 

1.7 

120 

CoFeB (5.4 nm, 2.5 nm) 

12 

-10 

87 

MgO (2.1 nm) 

* 

** 

4 

Ta (1.5 nm) 

0.75 

-5 

5.3 

Si 02 (500 nm) 

0 

- 

1 
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FIG. 6. A finite element simulation result showing the temperature profiles due to the Peltier effect without including the 
dissipation. Charge current of 150 /rA is sent through the MTJ resulting in a Peltier heating/cooling. If the temperature of 
one side of the junction; is anchored, the temperature change due to Peltier effect would be most prominent on the other side. 
Due to magneto-Peltier effect a different temperature at the top contact will be observed between P and AP configurations. In 
this figure, Sap= -99 /rV/K and 5'p= -145.2 /rV/K are used. To is taken as 290 K. 


The resistance of an MTJ is in general bias-dependent, i.e., the I — V curve is nonlinear. The I —V curves we obtain 
for device A in both the P and AP configurations are shown in Fig. 7. They are obtained by integrating the differential 
resistance dV/dl over current, which is acquired using a modulation technique where a small ac current (r.m.s 1 /iA) 
is superimposed on a dc bias in the lock-in measurements. For the P case, as can be seen from Fig. 7(a) (pink curve), 
the I —V characteristic is very close to a linear behavior, which is known as a special property for MgO-based MTJs. 
In this case the MTJ can be regarded as a normal resistor, where the Joule heating (power dissipation) at the junction 
(/ • V) is quadratically dependent on the current. For the AP case, however, the I — V curve is nonlinear especially 
at higher bias (grey curve), and the Joule heating effects would therefore deviate from a quadratic behavior, showing 
additional higher order dependences on the current. These additional heating effects thus appear as higher order 
thermoelectric signals at the thermocouple in the magneto-Peltier measurement. 

Suppose the I-V dependence can be expressed as: 

V{I) =a-I + h-I^ + c-I^ + d-I‘^ + ... (9) 

where a, 6, c and d are different order response coefficients. The signal at the thermocouple due to Joule heating can 
then be written as: 


Vthermocouple{I) — C ■ {a ■ + b ■ + C ■ + d ■ + .. .) (10) 

where C is a coefficient describing the efficiency of the conversion of heat generated in the MTJ to the voltage at the 
thermocouple, and is in a unit of A“^. Specifically, we are interested in the higher odd order responses shown up in 
the magneto-Peltier measurement, which come from the even parts of the I —V characteristics (related to coefficients 
b, d, ...), as can be directly seen from Eq. (9) and Eq. (10) . To discuss them separately, we extract the even parts 
of the I — V curves, by performing 

E(even) = ^[F(+/)+ !/(-/)] ( 11 ) 

on the I — V curves, and the results are shown in Fig. 7(b). Note that in Ref. [18], this part of signal is regarded as 
magneto-Seebeck signal as a result of the asymmetric heating induced temperature gradient across the MTJ. Here, 
however, using our finite element model we find this signal to be too large if explained in that way. We thus interpret 
the even parts as the intrinsic I — V characteristics of the MTJ, due to the inevitable difference between the two 
interfaces across MgO. 

For P configuration, the I —V curve follows an almost-linear behavior, and the even part is much smaller compared 
to the AP configuration. The higher harmonic signals generated from P configuration are all around 0 nV. From here, 
we only discuss the AP case. 

Now we calculate the odd Joule heating signal expected at thermocouple that is caused by this even part, by 
multiplying the coefficient C defined in Eq. (10). C can be determined in the much simpler parallel case, by comparing 
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FIG. 7. (a.) I — V characteristics of TMR for P (pink) and AP (grey) confignrations, respectively, (b) The even parts of the 
I — V characteristics extracted for these two configurations. 


the overall signal on thermocouple which is dominated by Joule heating with the I — V curve. Multiplying C and 
I with the grey curve in Fig. 7(b), we obtain an expected signal as the red solid line shown in Fig. S6. This is the 
expected signal on thermocouple that originates from the even I — V characteristics for the AP configuration due to 
dissipation. 

On the other hand, the odd order responses can be obtained by performing 

yodd order responses ^ 

= yi/ + y3/ _ y5f (_y7/) (^2) 

on the signals we obtained from lock-in measurements. The 1^7/ ig ignored here, as it is too small compared to other 
signals therefore omitting it does not change the overall analysis. Subtracting the linear response which we already 
know that originates from Peltier effect, we can get the odd non-linear signals, as plotted in Fig. 8 as black dots. It can 
be seen that the measured signals are very close to the expected values. It suggests that the nonlinear thermoelectric 
signals are indeed from Joule heating. 

In summary, the even component of the TMR I — V curve, that we clearly observe in our measurements, can 
produce Joule heating that shows higher odd order dependence with the current which mimics the magneto-Peltier 
effect. Note that this implies that Joule heating is different with opposite tunneling current direction, a similar 
observation as in Ref. [35], though the interpretation we have (based on careful harmonics analysis) is different 
from theirs (dissipation occurs at different contacts for opposite current direction). By comparing with our lock-in 
measurement signals detected on thermocouple we conclude that the higher odd harmonic signals that we observed 
are generated by this higher-order Joule heating effect. 



KmA) 

FIG. 8. Odd Joule heating signal that is expected to show on thermocouple in magneto-Peltier measurement configuration 
(red curve) and odd responses from magneto-Peltier measurement configuration that are calculated from lock-in measurements 
(black dots). Only the positive bias is shown. 








